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Abstract 

We present predictions for the J/ip and p meson production in the heavy ion ultraperipheral 
collisions (UPC) for the current energy 2.76 TeV at the LHC. Both total cross sections and 
cross sections with the neutron emission from one or both nuclei are presented. We also 
perform analysis of the RHIC /o-meson photoproduction data and emphasize importance 
of these data to test the current model for nucleus break up in UPC. 



1 Introduction 

Ultraperipheral collisions at the LHC provide a fine probe for studies of the small x dynam- 
ics in much cleaner environment than strong proton - proton (nucleus) collisions. They also 
provide an important bridge to the future studies at the e-A colliders. A detailed study of 
the UPC physics was performed by the UPC study group and reported in [1]. However, all 
these studies were performed for the planned energy of the heavy ion collisions at the LHC 
- y/s= 5.52 TeV. Hence, it is necessary to provide the predictions for the current energy 
y/s= 2.76 TeV in models, used in [T], as well as to explore long term advantages of having 
UPC data at two different LHC energies. Also, we receive requests from experimental 
groups to calculate the quantities which are more readily accessible experimentally - the 
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partial cross sections with dissociation of one or two nuclei with emission of a number of 
neutrons. 

The basic expressions for the cross section of production of vector mesons in ultrape- 
ripheral collisions are given in a number of papers, see e.g. [1]: 

da- AiA2^AiA2V 



dy 



n. 



■j/Ai iy)cr-yA2-^vA2 (y) + n^/A2 {-y)cr-yAi^vAi (-y) ■ (1) 



The quantity y = In^^ is the rapidity of the produced vector meson, and n^/Aiu) is the 
flux of photons with the energy = ^lQo emitted by one of nuclei (7_l is the Lorentz 
factor for colliding nuclei, and qo -is the photon momentum in the coordinate system of 
moving nucleus). 

The specific feature of the ultraperipheral collisions is rather large impact parameter 
b between interacting nuclei. Hence, the flux of photons, radiated by accelerated ion, can 
be with reasonable accuracy approximated by the equivalent photon spectrum of the point 
charge Z, moving with velocity (3 = pm/Ej^. To avoid the strong interaction of the colliding 
ions (which breaks coherence of the process), one has to restrict the impact parameters 
h of the interaction by choosing for example hmin ~ 2i?A or by suppressing the strong 
interactions by the factor Ps{h) = e:ic^[—aNNTAA{b)] with 

oo oo oo 

TAA{h)= [ dh' [ dzpA{b-b',z) [ pA{h',z)dz' 



describing overlap of colliding ions {pA{b,z) is the nuclear density normalized by the con- 
dition j d'^bdzp{b, z) = A). However, due to the large value of Z in UPC of heavy nuclei, 
there is a significant probability of the photon exchange with excitation of nuclei accom- 
panied by a subsequent decay with emission of neutrons [2] . When such additional photon 
exchange occurs in the final state of UPC with the vector meson photoproduction, this 
does not destroy the coherence of photoproduction, and sum over the final state of ions 
gives the total coherent photoproduction cross section. On the other hand, measurements 
by STAR at RHIC demonstrate that, using the Zero Degree Calorimeters (ZDC) which are 
installed in all detectors at RHIC and LHC, it is possible to select different channels, in 
particular, coherent photoproduction with excitation of ions in which neutrons are emitted 
by one or both nuclei. One can distingish several channels: Inln - emission of one neutron 
by each ion, XnXn - emission of several neutrons, Onln and OnXn excitation and decay 
only of one ion. To estimate these partial cross sections, we use the model suggested in [2] 
where it is shown that coherent photoproduction with additional electromagnetic nucleus 
excitation can be calculated modifying the flux of photons by impact parameter dependent 
factors P(j(h) which account for different channels (i=OnOn, Onln, Inln, OnXn, XnXn) 0. 
Hence, the final expression for the photon fiux used in our calculations is 



^Our thanks to S.Klein and J.Nystrand for the kindly providing their code for calculating the probabilities 
of electromagnetic excitation with subsequent neutron decay 
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7/a('^7 



vr 



dh- 



x) + 



KM 

7i 



where x 



and Kq,Ki are the modified Bessel functions. 



We provide predictions for the hght and heavy vector meson photoproduction with 
small transverse momenta of vector mesons. The dynamics of production of mesons, built 
of light quarks, and those, built of heavy quarks, is qualitatively different. Hence, we treat 
them separately. 



2 Production of light vector mesons 

Production of p-meson provides an important check of the modeling of the dynamics of 
UPC. Indeed, the cross section j + A^p + A is well understood theoretically. Though the 
Gribov space -time picture of high energy process of hadron - nucleus interaction differs 
qualitatively from low energy Glauber picture, numerically, the Gribov-Glauber theory 
which includes inelastic screening effects predicts very small (on the level of few percent) 
deviations from the Glauber formula for hadron-nucleus scattering. An additional element 
of the theory in the case of the photon projectile is presence of the non-diagonal transitions 
7 — )• /o' — )• /9. The generalized vector dominance model combined with Glauber approach 
describes well the data on p production at ~ lOGeV without free parameters j^. It 
also predicted [5] reasonably well absolute cross section of the p meson photoproduction for 
y/SNN =130 GeV as measured in the UPC at RHIC [6j. For high energy photoproduction 
off heavy nuclei correction to standard VDM+ Glauber model due to the account for the 
nondiagonal transitions is about 10-^ 15%. Since the uncertainty (combined experimental 
errors due to the restricted statistics and systematics) of the experimental cross sections 
measured by STAR exceeds this number almost by factor two, we use the standard Glauber 
model formula [7] for calculation of the cross section of coherent p meson photoproduction 
in the kinematics of UPC of heavy ions at RHIC and LHC : 



da.yN^pN{t = 0) , 

(^■yA^pA = / at 



dbe^i^-" / dzQ{b,z)e^\\ 



"-£f- f Q(b,z')dz' 2 



(3) 



Here q]_ = = tmin — t, —tmin = = is longitudinal momentum transfer in 

7-/9 transition. We also use Donnachie-Landshoff model [8] to calculate the total pN cross 
section and forward cross section of coherent p photoproduction off the nucleon target. This 
model describes reasonably well the data for 7p — pp at energies up to W^p k, 200 300 
GeV studied at HERA but can slightly overestimate the cross section at the very high 



^For the recent discussion see [5] 
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energies due to introducing the hard pomeron exchange with large intercept {an = 1-44) 
in such soft process as coherent photoproduction of p with small transverse momentum. 
Note, that we neglect the interesting effect of interference between production by left and 
right moving photons [9j since it gives a small correction in the cross section integrated over 
rapidities, pt- Our calculations for p meson photoproduction in UPC at RHIC and LHC 
are given in Tables [l]l3] and Figures [UlSl In Table [T] and Fig. [T] we compare cross sections 
obtained in the model described above with results presented by the STAR collaboration. 

It is worth noting that in the STAR experiment only the (XnXn) and (Inln) channels 
were measured. The p- meson were detected in the rapidity range — 1 < y < 1 and pt < 
150 MeV/c. The numbers reported by experiment for the total cross sections of various 
channels were calculated from the measured two cross sections using the the StarLight 
event generator to extrapolate to the larger \y\ and to calculate the ratio of the cross 
section in different channels. 

We find a reasonable agreement at lower energies but a significant discrepancy in all 
measured partial cross sections at ^/snn = 200 GeV. 
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f^lnln . mb 
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4.9 
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Table 1: Cross sections of the p meson photoproduction in gold-gold UPC at RHIC cal- 
culated in the Glauber model. Numbers in brackets present results obtained by STAR 
collaboration from analysis of the experimental data on the gold-gold UPC studies at 
RHIC. 

In particular, the data reported [10] do not show an increase of the cross section with 
increase of energy from = 130 GeV to \/sj^j^ = 200 GeV expected in all theoretical 

calculations. This is puzzling since the energy dependence of the cross section of p photo- 
production in UPC of heavy nuclei comes practically solely from the increase of the photon 
flux which is essentially model independent and, hence, very similar in all calculations. The 
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Figure 1: Comparison of the calculated cross sections for p photoproduction in the gold- 
gold UPC with the STAR experimental results. 
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Figure 2: Calculated rapidity distributions for different channels in p photoproduction in 
gold-gold UPC at RHIC. 
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Figure 3: Calculated rapidity distributions for different channels in p photoproduction in 
lead-lead UPC at LHC. 
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energy 


Partial cross section, mb 


a{XnXn, \y\ < 1) 


a{lnln, \y\ < 1) 


V^NN ~ ^"^^ 

GeV 


Glauber model 
STAR experiment 


18.6 
14.9 ±2.0 


1.87 
1.47 ±0.16 


V^iVAT = 200 
GeV 


Glauber model 
STAR experiment 


25.4 
14.5 ±2.0 


2.4 
1.07 ±0.16 



Table 2: Comparison of the Glauber model predictions with the results of measurements 
by STAR for p meson photoproduction in gold-gold UPC at \/s7V7V ~ ^00 GeV at RHIC. 
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5915 


8309 


^Lxn ' mb 
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1057 


'^XnXn ' mb 


261 


340 


<nln ' mb 


260 


306 


^Inln ' mb 


18.5 


21 



Table 3: Cross sections of the p meson photoproduction in PbPb UPC at LHC calculated 
in the Glauber model. 

cross section 7 + p— )-/>±pisa weak function of energy in the discussed energy interval 
and 7±j4— )-/9±^is expected to be even weaker function of the incident energy due to 
blackness of the interaction with heavy nuclei at small impact parameters. 

One can see from the Table [2] that the measured cross sections at two energies are 
practically equal while the calculated ones increase by a factor ~ 1.3, which coincides with 
an increase of the average photon flux. 

In Figl3]we present our predictions for the rapidity distributions of the p photoproduc- 
tion at energies ^fs^^ = 2.76 TeV and \fsj^j^ = 5.5 TeV at LHC. The total cross sections 
are given in Table [31 

Since da{'~fA — t- pA)/dt[t = 0) weakly depends on energy, combined studies at several 
energies and at different rapidities will provide a stringent test of the dynamics of the break 
up of nuclei due to e.m. interactions. 

Note in passing that our cross sections larger than the Starlight results since the 
Starlight MC code [11] uses expression for the total cross section of hadron - nucleus in- 
teraction given by classical mechanics which in the limit of large total VN cross sections 
and large A differs from the Gribov-Glauber model by a factor of about two. Recently one 
more update of predictions for total coherent cross sections in heavy ion UPC at RHIC 
and LHC has been published [12]. The calculations in [12] are based on use of Color Glass 
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Condensate (CGC) ideas and two versions of the color dipole model (IP-SAT [13j and 
IIM[T1]). Comparing to the data of STAR at \/sjvAr ~ ^00 GeV the authors found that 
IP-SAT model gives a larger cross section of p photoproduction but the result with IIM 
model appears to be close to that of the STAR experiment. So, cross section at the LHC 
energies predicted with this model are considered by the authors as preferable. However, 
as we show in this paper, the STAR measurements at y/sj^j^ = 200 GeV give the cross sec- 
tions practically equal to cross sections at ^/sJ^J^ = 130 GeV, and this doesn't agree with 
reasonable energy dependence, dictated by the increase of the photon flux with increase 
of 

V^NN- ^ result, the cross sections of p photoproduction in UPC at LHC, predicted 
in [12], appeared to be smaller than our. Note in passing that the color dipole models are 
usually used for description of the hard interactions. A rational to apply such a model for 
the description of the soft process of the p meson exclusive photoproduction is not clear, 
as opposed to the use of the (generalized) vector dominance model which is theoretically 
well justified in this case. 



3 J/i/j production 

The coherent production of J ftp in UPC is one of the direct ways of probing the pattern of 
the interaction of small color dipoles of average transverse size ~ 0.2-^0.25/m [15j. Due to 
their small size, the dominant mechanism of coupling of the cc to the nucleon is through a 
two gluon attachment which forms a start of the gluon ladder. In the dipole approximation 
the cross section of j + T ^ J/xjj + T process can be calculated in the Leading Order of 
PQCD [IlllIT! 



16r,.7r3 



dt 



as{p'^)xGT{x, p?' 



(4) 



Here Fee is width of the leptonic decay of J /if) and Gt{x,p'^) is density of gluons with 



fraction x 



of momentum of the target {s^t is the invariant energy for 7 — T 



scattering), and p'^ is the scale which we choose to be p'^ 



Applying EqH] to the 



proton and nucleus target, one can write the cross section for the photoproduction of J/ip 
off the nucleus in the form 



da^N^VNi^^, train) 

di 



M2 



Ga{-^,Q 



IGn{ 



^,Q^), 



2 imin 



dt 



dHdze''i^-^e-'i\\'p{h,z) 



(5) 



The key feature of the photoproduction of the hidden heavy flavor vector mesons off 
heavy nuclei is the gluon nuclear shadowing which is characterized by the ratio of the 
nuclear gluon density distribution G^(x,/i^) to the proton gluon distribution GAr(x,/i^). 



9 




10 





Figure 5: Calculated rapidity distributions for different channels in J /if) photoproduction 
in lead-lead UPC at LHC. 
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Figure 6: Calculated rapidity distributions for (OnOn) and (Onln) channels in J/^jj pho- 
toproduction in PbPb UPC at ^/sj^jq = 2.76 TeV. Shaded area shows uncertainty in 
accounting for the gluon nuclear shadowing. Left and right curves demonstrate one-side 
contributions. 
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Figure 7: Distributions on the transverse momentum transfer in coherent and incoherent 
J/V' photoproduction in PbPB UPC at i/sattv = 2-76 TeV. 
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The Ga/Gn ratio can be calculated within the theory of leading twist nuclear shad- 
owing [TU]( see [2D] for details of the model and references ) developed on the base of the 
Gribov's theory of inelastic shadowing, Collins's factorization theorem for hard diffraction 
and the DGLAP evolution equations. This model uses, as an input, nucleon diffractive 
parton distribution functions (pdf) which are available both in the NLO and in LO. In 
our study we select LO nucleon pdfs which give a reasonable description of the energy 
dependence of the elementary cross section. Also, we use the LO gluon density in pro- 
ton Gn{x,h = 2.5GeV^), recently found [18] by the Durham-PNPI group in the range 
10"^ <x < W-^ from the fit of cross section, given by (EqU]), to the HERA experimental 
data on 7 + p — J/tp + p. 

The results of calculations are presented in Figs. |l]-[6]and in the Table [H Since the 
estimate of shadowing is based on the parton distribution functions determined from the 
data, there is some uncertainty related to the experimental errors. In Fig. [6] we show how 
this uncertainty influence our estimates of the cross sections. 

Comparing our results with those of StarLight and [12], we find for J/-;/' photoproduction 
at LHC smaller cross sections. This is due to the accounting for the Leading Twist gluon 
shadowing which was not included in StarLight. The calculation of [12j included only 
higher twist shadowing in the color dipole model which is known to be much smaller than 
the leading twist shadowing in a wide range of x, see |22j for detailed comparison of the 
shadowing in the leading twist approximation and in the color dipole model. 

As we mentioned before, the coherent scattering is determined by the sum of lower and 
higher energy contributions. This makes it difficult to separate these contributions at any 
rapidity, except for y = where they are equal. One would like to determine the coherent 
cross section in as wide range of 'jA energies as possible, with higher energies allowing 
to probe the gluon density at lower x. For example, in PbPb UPC at -v/s^vAf ~ '^-'^^ 
TeV and rapidity of produced J/ip y ~ the corresponding x ~ 10~^, while, if one 
could extract contribution due to the higher energy photon to the the cross section of 
production J/^p with y 2 , it would be possible to extract gluon density at x w 1.5 • lO"'^. 
However, due to the rapid drop of the flux of photons with higher energy, this contribution 
is significantly lower than the cross section from the low energy photons (see Figl6|) which 
is determined by the gluon density at x « 10"^. It was argued [2] that measurement of 
the process of photoproduction with additional photon exchange enhances the contribution 
from the higher energy photons. We would like to emphasize that measuring at the same 
\fsj^iy two or more channels with J/ip at fixed rapidity, for example, (Onln) and (OnOn) 
or any other, one can easily separate high and low energy contributions, provided that 
the fluxes of photons modified by accounting for the neutron decay are calculated with 
reasonable accuracy. These fluxes and the procedure in whole can be tested in study of 
p photoproduction in PbPb UPC since the cross section of p photoproduction off nuclear 
target weakly depends on energy and can be calculated in the Glauber model. Also we 
would like to note that contribution of the incoherent photoproduction in the UPC at LHC 
which are usually accompanied by neutrons ^21j to the coherent cross section with detection 
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fOnOn, l^b 


10.1 


68 
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10.2 
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CXnXn, fJ-b 


3.7 


17 


548 
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fOnlnj fJ-b 


3.25 


16.8 


512 


739 


flnlnj l^b 


0.35 


1.6 


39 
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Table 4: Cross sections of Jip photoproduction in gold-gold UPC at RHIC and PbPb UPC 
at LHC calculated in the leading order LT shadowing model. 

of neutrons will be very small if one uses the low pt of produced J/tp (see Fig. [7|). 

If the measurements are done at two energies \/s= 2.76 TeV and -v/i= 5.52 TeV a si- 
multaneous analysis of two data sets will allow to extend model independent determination 
of nuclear ratio to x at least a factor of two smaller than in y = at -y/s= 5.52 TeV. 

4 Conclusions 

The studies of the J/ip production in the heavy ion UPC at the LHC will allow to measure 
coherent photoproduction cross section down to x w 10~^ using the channels without 
nucleus break up which dominate at the LHC energies. A cross check will be possible using 
channels with dissociation of nuclei. However, this would require sorting out situation with 
the /9-meson production at the RHIC energies. This requires both further experimental 
studies at RHIC and LHC as well as further theoretical studies. 

We would like to thank L.Frankfurt, V. Guzey, S.Klein, J.Nystrand, Yu.Gorbunov, for 
the discussions. This work was supported in part by the US DOE Contract Number DE- 
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of Science. 

References 

[1] A. J. Baltz, G. Baur, D. d'Enterria et al, Phys. Rept. 458, 1-171 (2008). 
[arXiv:0706.3356 [nucl-ex]]. 

[2] A. J. Baltz, S. R. Klein, J. Nystrand, Phys. Rev. Lett. 89, 012301 (2002). 

[3] M. Alvioh, C. Ciofi degh Atti, I. Marchino et al, Phys. Rev. C78, 031601 (2008). 

[4] L. Frankfurt, M. Strikman, M. Zhalov, Phys. Lett. B537, 51-61 (2002). 

[5] L. Frankfurt, M. Strikman, M. Zhalov, Phys. Rev. C67, 034901 (2003). 



15 



[6] C. Adler et al. [ STAR Collaboration ], Phys. Rev. Lett. 89, 272302 (2002). 
[nucl-ex7020600 4|. 

[7] T. H. Bauer, R. D. Spital, D. R. Yennie, F. M. Pipkin, Rev. Mod. Phys. 50, 261 
(1978). 

[8] A. Donnachie, P. V. Landshoff, Phys. Lett. B478, 146-150 (2000). |hep-ph/9912312l . 

[9] S. R. Klein, J. Nystrand, Phys. Rev. Lett. 84, 2330-2333 (2000). |hep-ph/9909"237l . 

[10] B. I. Abelev et al. [ STAR Collaboration ], Phys. Rev. C77, 034910 (2008). 
[arXi v:0712.332 [nucl-ex]]. 

[11] S. Klein, J. Nystrand, Phys. Rev. C60, 014903 (1999); 
|http: / /sourcefor ge.net /projects/upcstar light 

[12] V. P. Goncalves and M. V. T. Machado, Phys. Rev. C 84, 011902 (2011) 
[arXiv:1106. 3036 [hep-ph]]. 

[13] H. Kowalski, T. Lappi, R. Venugopalan, Phys. Rev. Lett. 100, 022303 (2008). 
[arX iv:0705.3047 [hep-ph]]. 

[14] E. lancu, K. Itakura, S. Munier, Phys. Lett. B590, 199-208 (2004). [hep- ph/0310338l . 

[15] L. Frankfurt, M. Strikman, M. Zhalov, Phys. Rev. Lett. 102, 232001 (2009). 
[arXiv:0811.0"368l [hep-ph]]. 

L. Frankfurt, V. Guzey, M. Strikman, M. Zhalov, JHEP 0308, 043 (2003). 

L. Frankfurt, L. Gerland, M. Strikman, M. Zhalov, Phys. Rev. C68, 044602 (2003). 

L. Frankfurt, M. Strikman, M. Zhalov, Phys. Lett. B540, 220-226 (2002). 

[16] M. G. Ryskin, Z. Phys. C57, 89-92 (1993). 

[17] S. J. Brodsky, L. Frankfurt, J. F. Gunion, A. H. Mueller and M. Strikman, Phys. Rev. 
D 50, 3134 (1994) |arXiv:hep-ph794 02283 1 . 

[18] A. D. Martin, C. Nockles, M. G. Ryskin and T. Teubner, Phys. Lett. B 662, 252 
(2008) (arXiv:0709.440lil [hep-ph]]. 

[19] L. Frankfurt and M. Strikman, Eur. Phys. J. A 5, 293 (1999) | arXiv:hep-ph/98 1 2322| . 

[20] L. Frankfurt, V. Guzey and M. Strikman, larXiv:1106. 20911 [hep-ph]. 

[21] M. Strikman, M. Tverskoy and M. Zhalov, Phys. Lett. B 626, 72 (2005) 

[22] L. Frankfurt, V. Guzey, M. McDermott, M. Strikman, JHEP 0202, 027 (2002). 



16 



